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Abstract 

The statistical coupling analysis of 768 p-glucosidases from the GHl family revealed 23 positions in which the amino acid 
frequencies are coupled. The roles of these covariant positions in terms of the properties of p-glucosidases were 
investigated by alanine-screening mutagenesis using the fall armyworm Spodoptera frugiperda p-glycosidase (SfPgly) as a 
model. The effects of the mutations on the SfPgly kinetic parameters {kcat/Km) for the hydrolysis of three different p- 
nitrophenyl p-glycosides and structural comparisons of several p-glucosidases showed that eleven covariant positions (54, 
98, 143, 188, 195, 196, 203, 398, 451, 452 and 460 in SfPgly numbering) form a layer surrounding the active site of the p- 
glucosidases, which modulates their catalytic activity and substrate specificity via direct contact with the active site residues. 
Moreover, the influence of the mutations on the transition temperature (7^) of SfPgly indicated that nine of the coupled 
positions (49, 62, 143, 188, 223, 278, 309, 452 and 460 in SfPgly numbering) are related to thermal stability. In addition to 
being preferentially occupied by prolines, structural comparisons indicated that these positions are concentrated at loop 
segments of the p-glucosidases. Therefore, due to these common biochemical and structural properties, these nine 
covariant positions, even without physical contacts among them, seem to jointly modulate the thermal stability of p- 
glucosidases. 
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Introduction 

In recent years, the search for new enzymes using improvements 
in sequencing technologies has resulted in a large collection of 
protein sequences. For example, the glycoside hydrolase family 1 
(GHl) groups have more than 5,000 P-glucosidase sequences in 
the CAZY database, only 274 of which have been marked as 
characterized to date [1]. Structural data - 41 crystaUographic 
structures are available - and biochemical characterization of GH 1 
P-glucosidases revealed these enzymes share the same fold, the (P/ 
Qt)8 barrel, and that their catalytic activity depends on a pair of 
glutamate residues, which act as acid/base and nucleophUe in a 
double substitution mechanism [1]. Moreover, their companions, 
R and Y residues, are involved in the modulation of the pK^ of the 
catalytic glutamates. These catalytic residues are highly conserved 
among P-glucosidases, except for the myrosinases [2,3]. Addition- 
ally, a network of hydrogen bonds formed by Q, H, W and E 
residues, which are placed and conserved in the active site of the P- 
glucosidases, modulates their substrate glycone specificity [4—7]. 
Finally, a set of variable residues, for which only the structural 
placement is conversed, forms the aglycone binding region for 
different P-giucosidases [8-10]. 

In parallel with the increase in sequence data, new methods to 
characterize the correlation between functions and structures of 
proteins have been developed that use special approaches to 



globally analyze protein sequences, and have revealed groups of 
residues that are joindy involved in determining functional 
properties. One of these methods is statistical coupling analysis 
(SCA), which, through covariation analysis of large multiple 
sequence alignments, is capable of identifying sets of residues that 
are important for protein folding [11], allostery [12], enzymatic 
activity and thermal stability [13]. Moreover, it was recendy 
demonstrated that sets of covariant residues, termed sectors, are 
important starting points for protein engineering [14]. 

Among the characterized GH 1 enzymes, the P-glucosidase from 
the fall armyworm Spodoptera frugiperda (SfPgly) has been extensively 
studied, including biochemical and site-directed mutagenesis of the 
active site residues involved in pH optimum modulation and 
substrate specificity [2,4,10,15,16]. Additionally, residues outside 
of the active site that affect SfPgly enzymatic activity via indirect 
contacts have been identified [1 7] . These findings make SfPgly an 
excellent model for analyzing the function of covariant residues of 
the GHl family. 

We applied SCA to an alignment containing 768 P-glucosidase 
sequences and identified 23 covariant positions. Using SfPgly as an 
experimental model, alanine residues were introduced at 18 of 
these covariant positions, and these single mutants were charac- 
terized for their thermal stability (TjJ and kinetic parameters {k^^J 
An,) for the hydrolysis of three different chromogenic substrates. 
Based on the results, a set of 1 1 covariant positions, which are 
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related to the enzymatic activity and form a layer surrounding the 
active site of the P-glucosidases, was identified. In addition, a 
second set of 9 covariant positions related to enzyme thermal 
stability, consisting of amino acid residues mosdy at the loop 
regions of this (P/ Q()jj barrel structure, was identified. 

Materials and Methods 

Identification of the covariant positions in |3-glucosidases 

SCA was performed as described previously [11,18] using a 
multiple sequence alignment containing 768 non-redundant P- 
glucosidases from the GH 1 family, which were retrieved from the 
PFAM server (http:/ /pfam.sanger.ac.uk). Site conservation (AG'"') 
and positional statistical coupling (AAG"°') parameters were 
calculated as previously described [19] using our own programs 
written in C/C^ [18]. Perturbations were considered significant 
when differences were present in at least 25% of the sequences in 
the alignment. A set of residue positions containing correlated 
conservations was obtained by clustering the AAG"'"' matrix using 
the Matlab (MathWorks) software package. 

Site-directed mutagenesis 

Site-directed mutants were constructed using the QuikChange 
site-directed mutagenesis kit (Stratagene, La JoUa, CA, USA) 
following the manufacturer's instructions. Wild-type Sfpgly cloned 
into the pAE vector [20] was used as a template for PGR reactions 
performed with mutagenic primer pairs, which are presented on 
the Table SI. The pAE vector coding for wild-type SfPgly, 
previously available, contains a T7 promoter, an ampicUlin 
resistance mark and the insert was cloned in the MM and Xhol 
sites. Mutation incorporation was verified by DNA sequencing. 

Expression and purification of recombinant and wild- 
type SfPgly 

NovaBlue (DE3) competent cells (EMD Millipore, Billerica, 
MA, USA) were transformed with pAE plasmids encoding wild- 
type or mutant SfPgly, plated on LB-agar containing ampicUlin 
(50 ng/mL) and grown at 37°G for 16 h. Single colonies were 
grown at 20°G in LB broth containing ampicUlin (50 [ig/ mL) until 
they reached an attenuance of 0.500 at 600 nm. Next, 0.4 mM 
(final concentration) isopropyl P-D- 1 -thiogsdactopyranoside 
(IPTG) was added for 16 h to induce recombinant protein 
expression, after which cells were harvested by centrifugation 
(4,000 x^, 20 min, 4°G) and frozen at -80°G. The pelleted cells 
were resuspended in 100 mM sodium phosphate pH 7.4 contain- 
ing 200 mM NaGl, 60 mM imidazole and 10% (v/v) glycerol and 
were lysed with 3 ultrasound pulses (15 s, output 10 using a 
Branson Sonifer 250 adapted with a microtip) with 1 min intervals 
in ice to avoid heating the samples. The supematants were 
recovered by centrifugation (13,200 g, 20 min, 4°G), and soluble 
recombinant proteins were incubated with 200 |J,L Ni-NTA 
Agarose (4°G, 1 h) (Qiagen, Hilden, Germany). The resin was 
pelleted and washed 5 times with 100 mM sodium citrate-sodium 
phosphate pH 6.0 containing 200 mM NaCl and 60 mM 
imidazole. Protein elution was performed using the same buffer 
containing 500 mM imidazole. Protein purity was verified by 
SDS-PAGE [21], and purified proteins were desalted using a 
minitrap G-25 (GE Healthcare, Upsala, Sweden) and stored at 
4°C. 

Protein concentration was determined by measuring the 
absorbance at 280 nm in 20 mM sodium phosphate pH 6.0 
containing 6 M guanidinium hydrochloride. Extinction coeffi- 
cients (8280 nm) wcrc Calculated based on the primary sequences of 
the wild-type and mutant Sfpgly proteins using the ProtParam 



server at ExPaSy (http://web.expasy.org/protparam/). The ex- 
tinction coefficients ranged from 2.149 to 2.053 (Table S2). 

Kinetic and thermal characterization of the mutant and 
wild-type SfPgly 

Enzyme kinetic parameters [kcax, and k^^t/K^ were 

determined for purified enzymes by measuring the initial rates 
(»()) of hydrolysis of at least 10 different concentrations of 
substrates, including /)-nitrophenyl-P-D-glucopyranoside, p-nitro- 
phenyl-P-D-fucopyranoside and /)-nitrophenyl-P-D-galactopyra- 
noside (Sigma, St. Louis, MO, USA) prepared in 100 mM sodium 
citrate - sodium phosphate buffer pH 6.0. Experiments were 
performed at 30°C. The hydrolysis of these substrates was detected 
following the formation of jft-nitrophenolate by absorbance at 
420 nm after addition of 250 mM sodium carbonate - sodium 
bicarbonate buffer pH 1 1 to the reaction samples. Kinetic 
parameters K„ and k^at were determined by fitting the »o and [S] 
data to the Michaelis-Menten equation using the Enzfitter 
software (Elsevier-Biosoft, Cambridge, UK). 

Differential scanning fluorimetry (DSF) experiments of both 
wild-t^pe and mutant SfPgly were performed using SYPRO® 
Orange solution (500-fold dilution) (Sigma, St. Louis, MO, USA). 
Melting studies were performed in optical tubes using a 7500 Real- 
Time System (Applied Biosystems, Foster Gity, GA, USA). The 
temperature gradient ranged from 25°G to 95°G with a slope of 
0.5% per step. The melting data were fitted according to recent 
literature [22], resulting in a theoretical Fitting processes were 
performed using EnzFitter software. 

Tridimensional modeling and computational structure 
comparison 

The tridimensional structure of SfPgly was homology-modeled 
based on the cr^'stallographic structure of the P-glucosidase from 
Neotermes koshunensis (PDB, 3VIK) using the Phyre2 software [23]. 
The sequences had 49% identity and 66% similarity. The 
structural models were visualized using the PyMOL molecular 
graphics system vl.l (Schrodinger, LEG). The distances between 
amino acid pairs were calculated using Deep View/ SwissPDB- 
Viewer v3.7 software [24]. The topology schemes of difiFerent P- 
glucosidases were manually drawn using information regarding 
the secondary structures visualized in PyMOL. 

Results and Discussion 

The SCA [19] using 768 sequences of GHl P-glucosidases 
revealed 23 covariant positions in their primary sequences (49, 54, 
57, 62, 98, 112, 143, 176, 188, 195, 196, 203, 223, 278, 309, 329, 
398, 445, 449, 451, 452, 456 and 460 in SfPgly numbering). 
Except for position 451, involved in the P-glucosidase's specificity 
for the substrate glycone [4—7], none had been experimentally 
correlated with any biochemical or biophysical properties of the P- 
glucosidases. 

Because SGA is based on the concept that covariant positions of 
a primary protein structure should be occupied by amino acid 
residues that are jointly involved in determining the same 
functional property [13], an experimental approach was designed 
to s(;arcli the identified set of covariant positions for sul)groups 
linked to the enzymatic activities or thermal stabilities of the GHl 
P-glucosidases. Therefore, SfPgly (GenBank code: AF 052729), a 
digestive enzyme from the fall armyworm Spodoptera fruff-perda that 
has been extensively studied [4,10,15], was chosen as a represen- 
tative P-glucosidase, and residues at its 23 covariant positions were 
separately replaced with alanine using a site-directed mutagenesis 
technique. Eighteen of these mutant enzymes were successfully 
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expressed as recombinant proteins in NovaBlue (DE3) bacteria 
and purified (Figure SI). The E451A mutant of Sffigly has been 
previously studied [4]. Mutant enzymes with replacements at 
positions 176, 329, 449 and 456 were not studied due to poor 
solubility. 

The enzyme kinetic parameters for the hydrolysis of three 
different substrates catalyzed by wild-type and the 18 mutant 
Sfpgly proteins were determined (Table 1 and Figure S2). 
Additionally, the transition temperatures [T^ for the denaturation 
of the mutant and wild-type SfPgly proteins were evaluated in 
Differential Scanning Fluorimetry (DSF) experiments (Table 2 and 
Figure S3). 

Mutations resulting in at least a 4-fold change in the tat/^m for 
the hydrolysis of at least two different substrates were considered 
replacements of residues relevant for the enzymatic activity 
because such variation corresponds to a AAG"'" higher than 3 kj/ 
mol, which is equivalent to the disruption or formation of one 
hydrogen bond [25] . Based on these criteria, only 4 mutations had 
no significant effects on the enzymatic activity, namely K49A, 
M57A, N112A and S445A (Table 3). 

Similarly, the mutational effect on the thermal stability of SfPgly 
was considered significant for Tj-^ changes higher than 2.5 K, 
given that in the DSF experiments, this variation range 
corresponds to four-fold changes in the ratio of native to denatured 
wild-type Singly. Based on that threshold, the mutations K49A, 
P62A, W143A, P188A, H223A, P278A, P309A, W452A and 
F460A were identified as replacements at positions involved in the 
thermal stability of SfPgly (Table 3). 



Table 3. Mutational effects 


on the catalytic activity (/Ccat/'^m 


ratio) and thermal stability (Aim) 


of SfPgly proteins. 




Mutation position 


kcaJKm ratio (mutant/wt) 




JT„ (K) 




NPpgIc 


NPpgal 


NPpfuc 




49 


0.5 


0.06 


0.7 


-7.2 


54 


0.02 


0.005 


0.02 


-0.2 


57 


0.77 


0.12 


0.8 


-1.0 


62 


0.09 


0.03 


0.18 


-5.5 


98 








-2.1 


112 


1.5 


0.19 


0.38 


-0.3 


143 


0.0004 


0.0003 


0.0005 


-7.2 


188 


0.0002 


0.001 1 


0.0012 


3.2 


195 


0.11 


0.01 


0.03 


0.1 


196 


0.06 


0.01 


0.06 


0.0 


203 


0.06 


0.04 


0.21 


0.8 


223 


0.30 


0.05 


0.075 


6.0 


278 


0.005 


0.02 


0.02 


5.8 


309 


0.15 


0.02 


0.02 


-7.1 


398 








nd 


445 


1.1 


1.5 


0.27 


-1.9 


452 








-2.8 


460 


0.0012 


0.0037 


0.0044 


-12.2 



Only mutational effects higher than a 4-fold change in the /CcatZ/^m ratio (0.25>/Ccat/^m ratio >4) were considered significant for the enzymatic activity. Mutational effects 
on the thermal stability were considered relevant only for ATm>2.5 K. -, no activity; nd, not determined; wt, wild-type. 
doi:l 0.1 371/journal.pone.0096627.t003 



Table 2. Transition temperatures (7^) for thermal 
denaturation of the wild-type and mutant SfPgly proteins. 



Enzyme 


7"m (K) 


wt 


319 


K49A 


331.6 


W54A 


318.8 


M57A 


317.9 


P62A 


313.5 


F98A 


316.9 


N112A 


319.3 


W143A 


311.8 


PI 88 A 


322.2 


G195A 


319 


Y196A 


319 


P203A 


319.8 


H223A 


325 


P278A 


324.7 


P309A 


311.9 


S445A 


317.1 


W452A 


316.2 


F460A 


306.8 



Standard deviations were lower than 0.5 K. 
doi:1 0.1 371 /iournal.pone.0096627.t002 
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IVFF 




Figure 1. Structural comparison of p-glucosidases sliowing tlie active site residues (red) and sector A positions (blue). Myrosinase 
from Sinapis alba (1 E6S); |3-glucosidase A from Paenibacillus polymyxa (1 EI4); p-glucosidase from Trichoderma reesei (3AHY); p-glucosidase Zmglu from 
Zea mays (1E56); p-glucosidase from Thermus thermophilus (1UG6); Human cytosolic p-glucosidase (2Z0X); SbDhr from Sorghum bicolor (1V03); p- 
glucosidase from Pyrococcus horikoshii (IVFF); p-glucosidase from Spodoptera frugiperda SfPgly. The distances between sector A and the active site 
residues are shorter than 4.5 A. The structures were visualized using PyiVlOL software. 
doi:1 0.1 371 /journal.pone.0096627.g001 



SC positions related to the catalytic activity in 
P-glucosidases 

Although no structural information was used in the SCA, visual 
inspection of the Sfpgly structural model revealed that residues at 
1 1 positions related to the enzymatic activity are directly 
connected to active site residues through covalent or non-covalent 
bonds. Therefore, those 1 1 residues form a single group indirecdy 
connected to each other through the active site residues and their 
non-covalent interactions with the substrate. This group was 
labeled sector A (Figure 1; Table 4). This observation is similar to 
previous SCAs of different protein families, which also showed 
groups of residues forming chains of interactions that connect 
distant points of their structures [19,26]. In the particular case of 
serine proteases, two sets of covariant positions were identified. 
The first set contained residues involved in the catalytic 
mechanism (including the catalytic triad), and the second set 
contained residues present in the SI pocket and was involved in 
substrate specificity [13]. 

Structural comparison showed that as in SfPgly, the placement 
of sector A positions in close contact to the active site is also 
observed for different fi-glucosidases from the GHl family 
(Figure 1; Table 4). In brief, residues at sector A positions contact 
two glutamate residues (essential for catalysis), a pair of residues 
(tyrosine and arginine) involved in the modulation of the pK.^ of the 
catalytic glutamate, and a group of residues that bind the 
substrates glycone and aglycone (Table 4). Thus, sector A residues 
might modulate several fi-glucosidase properties, from catalytic 
activity to substrate specificity by affecting the positioning and 
properties of their active site residues. Moreover, data presented 
here suggest that the active site of the fi-glucosidases is formed by a 



"layer" of highly conserved residues that interacts directly with the 
substrate and promotes its hydrolysis and is surrounded by a 
second "layer" formed by residues of the sector A positions 
(Table 4; Figure 1). 

Evidence that joint variation of residues at sector A positions is 
involved in substrate binding and specificity is found when 
comparing P-glucosidases and 6-phospho- fi-glucosidases, both 
groups belonging to the GHl family. Indeed, 99% of the 
sequences exhibiting Y at position 460 also present S at position 
451, and among these sequences, several were previously 
characterized as 6-phospho-P-glucosidases. Similarly, 99% of the 
sequences presenting E at position 45 1 present F at position 460, 
and several have fi-glucosidase activity. Thus, the identities of 
residues at positions 45 1 and 460 are linked and direcdy connected 
to enzyme specificity. In 6-phospho- fi-glucosidases, the replace- 
ment of E by S creates space in their active sites for the binding of 
substrates containing a phosphate group linked to the 6-OH of the 
glycone and also favors the formation of a hydrogen bond with this 
substrate. An additional hydrogen bond with the 6-phosphate 
group is formed by Y460 [31,32] (Figure S4). Therefore, changing 
the substrate specificity of these two groups of fi-glucosidases 
depends on replacements at both positions 45 1 and 460. A similar 
result was observed for the joint mutations C42A, C58A/V and 
S 1 95T at trypsin covariant residues, which converted this enzyme 
from a serine to threonine protease [13]. 

FinaUy, five residues, P62, G195, Y196, P278 and P309, 
(Table 3) were not included in sector A because no direct or even 
indirect interactions connecting them to the active site were 
identified. Thus, the replacements of those residues affected the 
Sffigly activity through currently unknown mechanisms. 
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1E4I 



1E56 



1\'03 
1E6S 




IVFF 



Figure 2. Distribution of sector S positions on tKie secondary structure of p-glucosidases. |3-glucosidase from Spodoptera frugiperda 
Sfpgly; p-glucosidase from Thermus thermophilus (1UG6); p-glucosidase A from Paenibacillus polymyxa (1E4I); p-glucosidase Zmglu from Zea mays 
(1E56); p-glucosidase SbDhr from Sorghum bicolor (1V03); myrosinase from Sinapis alba (1E6S); P-glucosidase from Trichoderma reesei (3AHY); Human 
cytosolic p-glucosidase (2Z0X); p-glucosidase from Pyrococcus horikoshii (IVFF). oc-Helices are represented by cylinders, p-strands by arrows and loops 
by lines. Sector S positions are shown as circles, whereas non-sector S positions are shown as stars. The symbols (circles or stars) in black indicate 
positions placed at loops, whereas white symbols mark positions at helices or strands. 
doi:l 0.1 371/journal.pone.0096627.g002 



SC positions related to the thermal stability of 
P-glucosidases 

The coupled positions in which replacements significantly afTect 
Sfpgly T„ (Table 3) are mostly occupied by prolines (4 out 9 
positions) or positively charged residues (K49 and H223). 
Conversely, only one proline residue was identified among the 
eight coupled positions not related to SfPgly thermal stability, and 
no charged residue is present among them. In addition, sequence 
comparison shows that positions 62, 188, 278 and 309 (Sfpgly 
numbering) are predominantly occupied by proline residues 
among family 1 P-glucosidase secjuences (frequencies higher than 
50%), whereas positively charged residues are dominant at 
positions 49 and 223 (90% and 63%, respectively). Moreover, 
structural comparison shows that positions related to thermal 
stability are concentrated at the loop segments of the P- 
glucosidases (8 out 9 for Sfpgly), whereas positions not related to 
thermal stabilitv' are mainly located in their ot-helices and P-strands 
(5 out 8 for SfPgly). Naturally, this distribution mainly results from 
the properties of the proline residues, which favor loop segments 
[33]. Indeed, these data are in agreement with previous 
observations demonstrating that proline and charged residues 
and loop length and mobilit}' are related to the thermal stability of 
proteins [34,35]. Thus, based on their similarities, coupled 
positions related to the thermal stability of P-glucosidases were 
labeled sector S (Figure 2). 

Therefore, in contrast to sectors observed in the serine protease 
[13] and sector A of P-glucosidases (discussed above), which are 
composed of residues contacting each other, sector S of the P- 
glucosidases is formed by positions with no physical interactions 



but is characterized by the prevalence of proline residues. Because 
of their placement in loops, sector S positions are jointiy related to 
the thermal stability of P-glucosidases. 

Supporting Information 

Figure SI SDS-PAGE of the purified wild-type and 
mutant SfPgly proteins. Mutant enzymes are identified by 

the residue and number of the covariant position followed by A, 
representing the replacement that they contain. Wt stands for wild 
type SfPgly. Details regarding the protein purification procedures 
are found in the Materials and Methods section. The gels (12% 
polyacrilamide) were stained using Coomassie blue R-250. 
Approximately 2 (xg of protein was loaded in each lane. MW, 
molecular weight markers. 
(TIF) 

Figure S2 The effect of substrate concentration on the 
activity of the wild-type and mutant SfPgly proteins. 

Mutant enzymes are identified by the residue and number of the 
covariant position in which A was introduced. The substrates are 
/(-nitrophenyl P-glucoside (NPPglc), /)-nitrophenyl P-galactoside 
(NPPgal) and /;-nitrophenyl P-fucoside (NPpfuc). Substrates were 
prepared in prepared in 100 mM sodium citrate - sodium 
phosphate buffer pH 6.0. Experiments were performed at 30°C. 
Dots are the experimental data. The continuous line represents the 
values calculated based on the best fitting of the experimental data 
into the Michaelis-Menten equation. Fitting was performed using 
EnzFitter software. 
(TIF) 
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Figure S3 Differential scanning fluorimetry of the wild- 
type and mutant SfPgly proteins. Samples (10 ng) were 
incubated with SYPRO stain, and the fluorescence was recorded 

at difierent temperatures. Mutant enzymes are identified by the 
residue and number of the covariant position in which A was 
introduced. Dots are the experimental data, whereas lines 
represent the calculated values produced based on the best fitting. 
Fitting process was performed using the software Enzfitter and the 
equation deduced on [22]. 
(TIF) 

Figure S4 Structural comparison between the active site 
of P-glucosidases (A) and 6-phospho P-glucosidases (B) 
from family GHl. A glucose unit (gray and red; A) is 
represented in the glycone subsite of the P-glucosidases SfPgly 
(green; Phyre2 model), Neokrmes koshunensis (blue; PDB 3VIF) and 
wheat Triticum aestivum (yellow; PDB 3AIQ), whereas a 6-phospho 
glucose unit (gray and red, B) is presented for the 6-phospho P- 
glucosidases from Streptococcus pneumonia (cyan, PDB 4IPN) and 
Lactococcus lactis (formely know as Streptococrux lactis) (majenta, PDB 
4PBG). Residues numbering was based on SfPgly. Dotted lines 
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